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We present results on (order v 2 ) QCD interquark potentials from SU(3) lattice 
simulations on volumes of up to 32 4 lattice sites at (3 — 6.2 and (3 = 6.0. Preliminary 
results on quarkonia spectra as obtained from these potentials are presented. 



1 Introduction 

Purely phenomenological or QCD inspired potential models have been proven 
to reproduce the observed charmonium (J/ip) and bottomonium (T) spectra 
remarkably well. The success of such models might be understood from 
QCD as a Schrodinger-Pauli type Hamiltonian with spin dependent (sd) and 
velocity dependent (vd) contributions — that can be parameterized in terms 
of seven independent scalar functions of the quark separation (potentials) - 
has been derived directly from the QCD Lagrangian [^, |3[]. By computing 
these potentials nonperturbatively on the lattice, the set of free parameters 
within the potential picture can be restricted. The present investigations 



1 Presented at the conference on quark confinement and the hadron spectrum (Confine- 
ment 96), Como. 
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correspond to a complete expansion of the QCD Lagrangian up to order 
1/m 2 in the heavy quark mass m or, alternatively, order v 2 in the relative 
velocity between the sources. The form of the Hamiltonian resembles that of 
NRQCD j| at this given order. The only additional assumption that enters 
is the instantaneous approximation, i.e. the potentials are assumed to have 
no explicit dependence on the interaction time. 

In our determination of the central potential as well as the relativistic 
corrections ||, we find short ranged Coulomb like contributions which had 
not been expected so far in two of the potentials. The main effect is that the 
J /if) "sees" a stronger effective Coulomb force than the T. 



2 The Hamiltonian 

Starting from a Foldy-Wouthuysen transformation of the Euclidean quark 
propagator in an external gauge field, the connection between the static 
interquark potential Vo(r) and Wilson loops can be derived. By perturbing 
the propagator in terms of the inverse quark masses mf 1 and rri^ 1 around its 
static solution, one arrives at the semi-relativistic Hamiltonian (in the CM 
system, i.e. p = Pi = — P2 and L = L x = L 2 ), 

H = Y,{mi + ^-- -f^] + V (r) + V sd (r, L, S 1; S 2 ) + V vd (r, p), (1) 

t=l V 2m i ) 

where the potential consists of a central part and sd and vd corrections [[l], 



1 m( mi J 2r 



with Rij = riTj/r 2 — <%/3 and, 

Kd(r,p) = U— 2 + ^)(V 2 V (r) + V 2 V a (r)) (3) 
8 \mf 171% J v ' 



1 , „ . A 1 



{PiiPj-i Sij\w ^ X/ ™2 {PiiPji Tij} w 



2 



with Sij = SijV b (r) - RijV c (r) and T {j = 5ijV d (r) - R i:j V e {r). The symbol 
{•, •, -}w = H - ) {•) •}} denotes Weyl ordering of the three arguments. V[ - 
V4 are related to spin-orbit and spin-spin interactions, Vb-V e to orbit-orbit 
interactions and the Darwin-like term that incorporates V 2 V a modifies the 
central potential. V[ - V4 and V 2 V a - V e can be computed from lattice 
correlation functions of Wilson loop like operators in Euclidean time. Pairs 
of the potentials are related by Lorentz invariance to the central potential 
(y 2 ' _ V ; = V h - 2V C = \V Q - $V{, V c + 2V e (r) = -§V£) g, §, such that 
only six out of them are truly independent. 



3 Results 

We find the sd potentials V^', V 3 , V4 to be short ranged and to agree with 
tree level lattice perturbation theory. The agreement is even better when a 
running coupling in momentum space is allowed. Apart from a constant part, 
which accounts for the string tension, the spin-orbit potential V[ contains an 
attractive Coulomb-like piece which is not expected from perturbation theory 
and makes up about 20 % of the Coulomb part of the central potential. To 
leading order minimal surface approximation, V a does not contain a long 
range component while the short range part is expected to vanish from weak 
coupling perturbation theory. We, however, find it to behave like — 1/r which 
means that the Coulomb force within cc systems is enhanced by about 20 % 
in comparison to the bb case. V\, - V e are found to agree qualitatively with 
the combined long and short distance expectations. By neglecting the very 
short range behaviour (r < 0.15 fm), which is affected by running coupling 
effects, the data can be parameterized as follows: 

e h d 

V = KT-- V 2 K = — , V{ = -R - — 

Vi = - 2 , V 3 = 3 -^, V, = 87ra5 3 (r), (4) 

2 e KT e KT KT KT 

with a = e — d. V , V a , V& and contain unphysical self energy constants 
which have to be subtracted. The parameter values are: e = 0.321(7), d = 
0.065(11) and h = 3.75(31)«. 
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Figure 1: The bottomonium spectrum (in GeV). Horizontal lines are experimental 
results. 

4 Spectroscopy 

Two free QCD parameters have to be determined from experiment, i.e. from 
a fit to the T and J /if) spectra, namely the heavy quark mass, m& or m c , and 
a scale which sets the initial condition for the running of the QCD coupling, 
e.g. the parameter Since we have not included sea quark effects into 

2 Since the central potential is unbounded, the zero point energy is in principle arbitrary. 
However, it turns out that a zero point energy C in Vq will give rise to counter terms in 
Vb and Vd which can be absorbed into a redefinition of the quark mass J3|: m — ► m + C/2. 
We choose C = by demanding Vd(0) — 0. Apart from this, the constituent quark masses 
of Eq. ^] can differ from the kinetic masses appearing in the dispersion relation to this 
order in 1/m. We find that allowing for such a difference does not substantially improve 
the agreement with experiment. 
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our lattice simulation yet, for the time being, we allow for a free Coulomb 
coefficient e. The optimal values of the fit parameters turn out to be e = 
0.495, ■ s /k « 430 MeV, m b « 4.79 GeV and m c m 1.41 GeV with an error 
of about 0.1 on e. Simulations of full QCD || suggest that values e ~ 0.40 
are realistic in a world with three light sea quarks. In Fig. [I] the bb spectrum 
from our best fit (e = 0.495) is compared to experiment (horizontal lines) 
as well as the result with e = 0.4. Fits to improved parametrizations that 
incorporate the QCD running coupling are in preparation. 
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